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In this paper, we construct a new class of charged rotating dilaton black brane solutions, with 
complete set of rotation parameters, which is coupled to a nonlinear Maxwell field. The Lagrangian 
of the matter field has the form of the power-law Maxwell field. We study the causal structure 
of the spacetime and its physical properties in ample details. We also compute thermodynamic 
and conserved quantities of the spacetime such as the temperature, entropy, mass, charge, and 
angular momentum. We find a Smarr-formula for the mass and verify the validity of the first law 
of thermodynamics on the black brane horizon. Finally, we investigate the thermal stability of 
solutions in both canonical and grand-canonical ensembles and disclose the effects of dilaton field 
and nonlinearity of Maxwell field on the thermal stability of the solutions. We find that for a < 1 , 
charged rotating black brane solutions are thermally stable independent of the values of the other 
parameters. For a > 1 , the solutions can encounter an unstable phase depending on the metric 
parameters. 

PACS numbers: 04.70.Bw, 04.30.-w, 04.70.Dy 


I. INTRODUCTION 

Historically, the presence of a scalar field in the context of general relativity dates back to Kaluza-Klein theory [l| . 
Kaluza-Klein theory was presented in order to give a unified theory describing gravity and electromagnetic. It may 
be considered as a pioneering theory to the string theory. Afterwards, in Brans-Dicke (BD) theory, the notion was 
seriously taken into account Q ■ The root of this theory can be traced back to Mach’s principle encoded by a varying 
gravitational constant in BD theory. The variety of gravitational constant is indicated in this theory by considering a 
scalar field non-minimally coupled to gravity. BD theory is known as a non-quantum scalar-tensor theory 13]. From 
quantum viewpoint, the scalar field called dilaton field is emerged in the low energy limit of string theory [J]. 

On the other side, there are strong evidences on the observational sides that our Universe is currently experiencing a 
phase of accelerating expansion [5] . This acceleration phase cannot be explained through standard model of cosmology 
which is based on the Einstein theory of general relativity. One way for explanation of such an acceleration is to add 
a new unknown component of energy, usually called “dark energy”. Another way for challenging with this problem 
is to modify the Einstein theory of gravity. In this regards, some physicists have convinced that Einstein theory of 
gravitation cannot give a complete description of what really occurs in our Universe and one needs an alternative 
theory. As we mentioned above, string theory, in its low energy limit, suggests a scalar field which is nonminimally 
coupled to gravity and other fields called dilaton field [ 5 ] . Consequently, dilaton gravity as one of the alternatives of 
Einstein theory have encountered high interest in recent years. String theory also suggests dimensions of spacetime to 
be higher than four dimensions. It seemed for a while that dimensions higher than four should be of Planck scale, but 
recent theories expresses that our three-dimensional brane can be embedded in a relatively large higher dimensional 
bulk which is still unobservable j§j, In such a scenario, all gravitational objects including black objects are higher 
dimensional. The action of dilaton gravity usually includes one or more Liouville-type potentials which can be justified 
as trace of supersymmetry breaking of spacetime in ten dimensions. 

Many authors have explored exact charged black solutions in the context of Einstein-dilaton gravity. For instance, 
asymptotically flat black hole solutions in the absence of the dilaton potential are studied in Refs. fa-tlll]. In recent 
years, AdS/CFT correspondence have made asymptotically non-flat solutions more valuable. On the other hand, 
study of asymptotically non-flat and non-AdS solutions extends the validity range of tools already applied and tested 
in the cases of asymptotically flat or AdS spacetimes. Many efforts have been made to find and study static and 
rotating nonflat exact solutions in the literature. Static asymptotically nonflat and non AdS solutions and their 
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thermodynamics were explored in [IMzl- Different aspects of exact rotating nonflat solutions have also been studied 

in mill. 

In the present paper we extend the study on the dilaton gravity to include the power-law Maxwell field term in the 
action. The motivation for this study comes from the fact that, as in the case of scalar field which has been shown 
that particular power of the massless Klein-Gordon Lagrangian shows conformal invariance in arbitrary dimensions 
[28 ], one can have a conformally electrodynamic Lagrangian in higher dimensions. It is worth mentioning that 
Maxwell Lagrangian, F llv F lllJ is conformally invariant only in four dimensions, while it was shown that the Lagrangian 
(F^Fn (n+1)/4 is conformally invariant in (n + l)-dimensions }28j j. In other words, this Lagrangian is invariant under 
the conformal transformation —> G 2 ^,, and A^ —» A^. The studies on the black object solutions coupled to a 

conformally invariant Maxwell field have got a lot of attentions in the past decades [29(-[33j. Thermodynamics of 
higher dimensional Ricci flat rota ting black branes with a conformally invariant power-Maxwell source in the absence 
of a dilaton field were studied in (32] ■ Recently, we explored exact topological black hole solutions in the presence of 
nonlinear power-law Maxwell source as well as dilaton field 33j. Since these solutions (|HJ) are static, it is worthwhile 


to construct the rotating version of the solutions. Charged rotating dilaton black branes in dilaton gravity and their 
thermodynamics in the presence of linear Maxwell field and nonlinear Born-Infeld electrodynamics have been studied 
in [27] and [HI respectively. Till now, charged rotating dilaton black objects coupled to a power-law Maxwell field has 
not been constructed. In this paper, we intend to construct a new class of (n + l)-dimensional rotating dilaton black 
branes in the presence of power-law Maxwell field, where we relax the conformally invariant issue for generality. Of 
course, the solution do exist for the case of conformally invariant source. We find that the solutions exist provided one 
takes Liouville-type potentials with two terms. In the limiting case of Maxwell field, one of the Liouville potentials 
vanishes. We shall investigate thermodynamics as well as thermal stability of our rotating black brane solutions and 
explore the effects of the dilaton and power-law Maxwell fields on thermodynamics and thermal stability of these 
black branes. 

The outline of this paper is as follows. In the next section, we present the basic field equations of Einstein-dilaton 
gravity with power-law Maxwell field and find a new class of charged rotating black brane solutions of this theory and 
investigate their properties. In section uni we obtain the conserved and thermodynamic quantities of the solutions 
and verify the validity of the first law of black hole thermodynamics. In Sec. IIV1 we study thermal stability of the 
solutions in both canonical and grand-canonical ensembles. We finish our paper with closing remarks in the last 
section. 


II. FIELD EQUATIONS AND ROTATING SOLUTIONS 


We consider an (n + l)-dimensional (n > 3) action of Einstein gravity which is coupled to dilaton and power-law 
Maxwell field, 


where 1Z is the Ricci scalar, $ is the dilaton field and U(<I>) is the dilaton potential. Here F = F^F^ where 
Ffiu = d[p,A v ] is the electromagnetic field tensor and A M is the electromagnetic potential, while p and a are two 
constants that determine nonlinearity of electromagnetic field and coupling strength of the dilaton and electromagnetic 
fields, respectively. The well-known Einstein-Maxwell dilaton theory corresponds to the case p = 1. The last term 
in (HD is Gibbons-Hawking boundary term which is added to the action in order to make variational principle well- 
defined. We show the metric of manifold A4 with g^, with covariant derivative V^. The metric of the boundary 8A4 
is "fab and 0 = 7 o &0 ab is the trace of the extrinsic curvature of the boundary 0 ab . By varying action CD with respect 
to the gravitational field g the dilaton field $ and the gauge field A one can obtain equations of motion as 


V 2 $ — n ~ ^ ^ -Aup3>/(n-l) 1 _ r\p 

8 5 $ + 2 1 ’ ’ 

(V 4ap$/(n “ 1) (-.7 r ) p “ 1 F'“'') = 0. 


4 

n — 1 


+ 2pe- 4ap * / ( n - 1 \-F) p - 1 F^ x F u \ 


( 2 ) 

( 3 ) 

( 4 ) 


Here we seek for higher dimensional rotating solutions of field equations (|2 ]l -(|4 |) . We know that rotation group in 
(n+ l)-dimensions is SO(n). Thus the number of independent rotation parameters for a localized object is k = [n/2] 
where [a:] is the integer part of x. Note that k is equal to number of Casimir operators. Therefore, the metric of 
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(n + l)-dimensional rotating solution with cylindrical or toroidal horizons and complete set of rotation parameters 
can be written as 35] 


ds 2 = -f(r)(Edt-^2 a idfa j + j^R 2 (r)(CLjdt - El 2 dfa ) 


i= 1 


r 2 k dr 2 r 2 

~ ajdfaf + —rr + -pR 2 (r)dX 2 

i<j J ^ > 


k 2 


= i+ £! 


i =1 


(5) 


where cqs are fc rotation parameters and l has the dimension of length which is related to the cosmological constant 
A for the case of Liouville-type potential with constant <h. The angular coordinates are in the range 0 < fa < 2n 
and dX 2 is the Euclidean metric on the (n — k — l)-dimensional submanifold with volume E ra _fc_i. Integrating the 
Maxwell equation 0 leads to 


Ft r — 
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-pS 2 = 0 , 


( 6 ) 


where q is an integration constant related to the electric charge of the brane. Substituting (0 and 0 in the field 
equations m and 0, we find the following differential equations 


(7) 


r _ (n-l)fR' _ (n — 1) f _ 4/£^ _ 2(n-l)fR! _ (n - 1) fR" _ E(<1>) 
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4q p<3? 

(n-3) //' i?7'\ /i?" 2(n-l)/i? / 2 P pq 2 P e ^-^p-^ 

~ V7 + 7T J ” 7? 7? 


(13) 


where the prime denotes derivative with respect to r. In order to construct exact rotating solutions of the theory 
given by action 0, the functions f(r), R(r) and <E>(r) should be determined so that the system of equations (I71)- lfl3l) 
are satisfied. To do that, we make the ansatz [27j 


R(r)=e 2a ^ r)l{n ~ 1 \ (14) 

and take the potential of the Liouville-type with two terms, namely 

y($) = 2A 1 e 2Cl<I> + 2Ae 2C2<I> , (15) 

where Ai, A, (i and C ,2 are constants. Using m and m , one can easily show that equations m-m have solutions 
of the form 


f(r) = 


2Ab 2 ' Y (l + a 2 ) 2 r 2 ( 1-7 ^ m 2 p p(l + a 2 ) (2p - 1) q 2p 

(n - 1) (n - a 2 ) r («-i)(i-7)-i + ’ 


(16) 


<l>(r) 


(n-l)a 

2(1 +a 2 ) n 



(17) 


where b and m are arbitrary constants, 7 = a 2 /(a 2 + l), T = (n—2p+a 2 )/(2p— l)(l+a 2 ) and II = a 2 + (n — 1 — a 2 ) p. 
The above solutions will fully satisfy the system of equations provided, 


2p (n — 1 + a 2 ) 2a 

{n — 1 ) (2p — 1 ) a ’ 2 n — 1 


Ar 


2 P_1 ( 2 p — 1 ) (p — 1 ) a 2 q 2p 


Ub 


2(n-l )p 
2 p-1 


It is worthwhile to note that the necessity of existence of the term 2Aie 2 ^ 1$ in scalar field potential is due to 
nonlinearity of Maxwell field and this term vanishes for the case of linear Maxwell field where p = 1 (27} ■ Note that 
in our solutions A remains as a free parameter which plays the role of the cosmological constant. Another thing to 
notice is that although these rotating solutions are locally the same as those found in [33] with flat horizon {k = 0), 
they are not the same globally. One may also note that in the particular case p = 1 these solutions reduce to the 
[n + l)-dimensional charged rotating dilaton black branes presented in (23] • The parameter m in Eq. (fTBl) is the 
integration constant which is known as the geometrical mass and can be written in term of the horizon radius as 


m(r+) 


2 p p (2p — 1) q 2p b 21 nr+ 7)( ” +1) 1 

(l- 7 )Tnfe 2 ('*- 2 )TP/( 2 i’- 1 VJ + / 2 (l-7 ) 2 (n-a 2 )’ 


(18) 


where r+ is the positive real root of /(r+) = 0. One can easily show that the vector potential A^ corresponding to 
the electromagnetic tensor © can be written as 




(2p+l-n) 7 

qb ( 2 ,- 1 ) 

t7 


(H-art) 


(no sum on j). 


(19) 


Here we pause to give some remarks about the value of A M at infinity that apply some restrictions on the value of p 
and a. This discussion is given here because it is necessary to know these restrictions for next studies. In section (ED, 
we will show that the mass of solutions is dependent on the value of A M at infinity. Thus, the value of electromagnetic 
vector potential should be finite at infinity so that we have finite mass. In order to guarantee this behavior T should 
be positive i.e. 


n — 2p+a 2 
{2p — l)(l+a 2 ) 


The above equation leads to the following restriction on the range of p, 
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n + a 2 
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On the other hand, the effect of m in metric function should vanish in spacial infinity. This fact leads to a restriction 
on a 

a 2 < n - 2. (22) 

Therefore, one can sum up the constraints m and (EH) as follows: 

1 n o . 

for - < p < — , 0 < a 2 < n — 2, (23) 

Zj Z 

Tl 

for — < p < n — 1, 2p — n < a 2 < n — 2. (24) 

It is worth mentioning that the solution is always well-defined in the above allowed ranges of a and p. In continue of 
this section we discuss properties and asymptotic behaviors of the solutions in allowed ranges of p and a. 


A. Asymptotic behaviors and properties of the solutions 

Now, by considering the restrictions (E?l) and (EH) on p and a, we are ready to discuss the behavior of our solutions 
both in the vicinity of r = 0 and infinity. First, one should note that in allowed ranges of p and a, the charge term in 
metric function disappears at infinity as the mass term does. This behavior is also seen in the special cases of p = 1 
or a = 0. Therefore, the first term in f(r) 1 (1161) 1 determines the behavior of it at infinity. Since we are interested in 
the solutions that go to infinity as r — > oo, we assume A < 0 and take it in the standard form A = — n(n — l)/2 1 2 . It 
is also notable to mention that the spacetime is neither asymptotically flat nor AdS due to existence of dilaton field. 
About r = 0, the term including q in metric function is dominant. This term is always positive in permitted ranges 
of p and a because II, T > 0 in these ranges. Therefore as in the case of Reissner-Nordstrom black holes, we have 
timelike singularity and there are no Schwarzschild-type black hole solutions. 

We continue our discussions about properties of the solutions by looking for curvature singularities. Since essential 
singularities are located at divergencies of Kretschmann scalar, we seek for these for our solutions. It is easy to 
show that the Kretschmann scalar R ll v\nR fJ,uXl ‘ i diverges at r = 0 while it is finite for r / 0 and goes to zero as 
r — > oo and therefore there is an essential singularity located at r = 0. Although the location of event horizon cannot 
be determined analytically by using f(r), fortunately we can get more insight about the solutions by calculating 
temperature corresponding to event horizon. The temperature and angular velocity of the horizon can be obtained 
by analytic continuation of the metric. The analytical continuation of the Lorentzian metric by t —>■ it and a —> ia 
yields the Euclidean section. In order that Euclidean metric is regular at r = r+, one should identify r ~ r + /3 + and 
4>i ~ (fri + where /3+ and fl^s are the inverse Hawking temperature and the ith component of angular velocity 

of the horizon. Then, temperature and ith component of angular velocity can be computed as 


f'{r+) (1 + a 2 ) ( nb 2l r\ 2l 2 p p (2p — 1) q 2p 1 

47rS 47 tS | l 2 JJ52(n-2)7p/(2p-l) r J+( 1 -7)(™-l) J ’ 


(25) 

(26) 


Numerical calculations show that temperature vanishes at event horizon for extreme black brane solutions. Therefore, 
one can see from Eq. (1251) that we have extreme black brane if 


TTT-ext 


( a 2 + 1 )nb 2j 


(1 + Q: )T 
(n-a 2 ) \ rext 


1 


(27) 


or 


„ U n b2l(np-l)/(2p-l) 

2p _ _ T+(l- 7 )(n+l)-l 

yext 2 Pp (2 p - 1) l 2 ext 


(28) 


The solutions have two inner and outer horizons located at r_ and r + , provided the charge parameter q is lower than 
q ex t or m is greater than ?Ti ex t an d a naked singularity if q > q ex t or m < m ex t (see Fig. [lj. Note that there is a 
relation between m ext and q exi as 
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/ (a 2 + l)n& 27 ((n - a 2 ) + (1 + a 2 )T) \ / 2 p p(2p-l)l 2 2p 

mext “ ( Tl 2 ( n-a 2 ) J \nnb 2 7( n P-V/( 2 P-V 9ext 

which reduces to extremal mass obtained in [27] for linear Maxwell field (p = 1) and to one obtained in [43[ in the 
absence of dilaton field a =7 = 0 and with linear Maxwell field. Finally, it is noticeable to mention that there is 
a Killing horizon in addition to event horizon for rotating solutions in our Einstein-dilaton gravity as in the case of 
rotating black solutions of the Einstein gravity. It is easy to show that the Killing vector field 


(l- 7 )(n + l)-l 


(29) 


X = Cx, 

k 

X = + ( 3 °) 

2=1 

is the null generator of the event horizon, where k denote the number of rotation parameters [36j and C is a constant 
that we will fix it in next section. The Killing horizon is a null surface whose null generators are tangent to a Killing 
field. 



FIG. 1: The behavior of f(r) versus r with l = b = 1, q = 0.5, a = y/2, n = 5, p = 2, r ex t = 1. In this case m ex t = 60.00 and 

<?ex t = 1.06. 


III. THERMODYNAMICS OF BLACK BRANES 


In this section, we discuss thermodynamics of rotating black brane solutions in the presence of power-law Maxwell 
field. Since the discussion of thermodynamics of these solutions depends on the calculation of the mass and other 
conserved charges of the spacetime, we first compute these conserved quantities. The way we use for calculating 
conserved quantities is counterterm method. This method is a well-known method for Asymptotically Ads solutions 
to avoid divergencies in calculation of conserved quantities which is inspired by AdS/CFT correspondence [37:]. This 
method can also be used for dilaton gravity and in the presence of Liouvill-type potential where the spacetime is not 
asymptotically AdS [22|, [27j, (13■ Since boundary curvature of our spacetime is zero (R a bcd( 7 ) = 0), the counterterm 
for the stress energy tensor should be proportional to j ab . We find the finite stress-energy tensor in (n-fl)-dimensional 
Einstein-dilaton gravity with Liouville-type potential as 


T 


ab 


1 

8-7T 



0 7 a6 


n — 1 

^eff 



(31) 


where Z e ff is given by 


, 2 = (n- l)(a 2 -n) 
eff V($) 


(32) 
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The first two terms in m are the variation of the action (HD with respect to 7 ab , and the last term is counterterm 
which removes the divergences. Note that in the absence of the dilaton field (a = 0), we have K(<1>) = 2A, and the 
effective l 2 s of Eq. (1321) reduces to l 2 = — n(n — 1)/2A of the AdS spacetimes. In order to compute the conserved 
charges of the spacetime, one should first choose a spacelike surface B in dAA with metric 07 ,-, and write the boundary 
metric in ADM (Arnowitt-Deser-Misner) form 

'y ab dx a dx a = — N 2 dt 2 + &ij (dip 1 + V l dt ) (dpj + V^dt) , 

where the coordinates p l are the angular variables parameterizing the hypersurface of constant r around the origin, 
and N and V 1 are the lapse and shift functions respectively. Then, the quasilocal conserved quantities associated 
with the stress tensors of Eq. ED can be written as 


Q(0 = [ d n -\^T ab n a i b , (33) 

Jb 

where <7 is the determinant of the metric <t,j, £ and n a are the Killing vector field and the unit normal vector on 
the boundary B. In order to calculate quasilocal mass and angular momentum, one should choose boundaries with 
timelike (£ = d/dt) and rotational (? = d/dp ) Killing vector fields i.e. 


and 


M = 




(34) 


J = 


d n ~ x p^T ab n\ b , 


(35) 


provided the surface B contains the orbits of <7 (1551) and (1551) are conserved mass and angular momenta of the black 
hole surrounded by the boundary B. It is remarkable to mention that although mass and angular momenta do not 
depend on specific choice of foliation B within the hypersurface dA4 , they are dependent on location of boundary B in 
the spacetime. Taking into account the cylindrical symmetry of the rotating black brane with k rotation parameters 
along the angular coordinates 0 < (f>i < 2i r, we denote the volume of the hypersurface boundary at constant t and r 
by V n -\. Then, the mass and angular momentum per unit volume V n -\ of the black branes can be calculated through 
the use of Eqs. (l34| and (1551) . We find 


{,(«— 1)7 
167rZ n_2 



a 2 )z? + a 2 — 1 
1 + a 2 


TO, 


(36) 


J)(n- 1)7 
16nl n - 2 



(37) 


As one can see from ED, the angular momentum per unit volume is proportional to a,iS and therefore it vanishes 
if di = 0 (S = 1). Thus, it is physically reasonable to consider a^s as rotational parameters of the spacetime. The 
last conserved quantity of our solutions is electric charge. Electric charge can be obtained by calculating the flux of 
the electric field at infinity. By projecting the electromagnetic field tensor on specific hypersurfaces, we can find the 
electric field as E M = g^Pe~ iap ^^ n ^ 1 ' ) (—F) p_1 F pv u v where u v is normal to such hypersurfaces. The components of 
u v are 


u 


0 


1 

N’ 


= 0 , 


V/_ 

~N’ 


(38) 


where N and V 1 are the lapse function and shift vector. Eventually, the electric charge per unit volume I 4.-1 can be 
calculated as 


Q = ——— 

^ 47W"- 2 ’ 

q=2 p - 1 q 2p ~ 1 . (39) 

One may note that q = q for p = 1 [27j . 

Now, we calculate thermodynamic quantities entropy S and electric potential U. Entropy of almost all black 
solutions including ones in Einstein gravity typically obeys the so called area law [39j, |4Ci} . Dilaton black solutions are 







not exceptions (see for instance [13, HU). Thus, we can calculate the entropy per unit volume V n -i of our rotating 
black brane as 


S = 


4 l n ~ 2 


(40) 


The electric potential U, can be calculated through the definition i filj 

u = A„r\r^oo - A^r \ r=r+ , (4i) 

where \ is the null generators of the event horizon given by Eq. ®. Therefore, the electric potential may be obtained 
as 


U = 


Cqb 


(2p+l-n) 7 

(2p-l) 


STr T 


(42) 


Here, we are ready to seek for satisfaction of thermodynamics first law. First, we should obtain the mass M in terms 
of extensive quantities S, Q and J. Using Eqs. (I36l) - (l39l) and the fact that f(r+) = 0, we receive 


M(S,Q, J) 


((n — a 2 )Z + a 2 — l) J 
(n — u 2 )l\/Z(Z — 1) ’ 


(43) 


where J =\/J2i J? and % 


c? which is the positive real root of the following equation: 


( a 2 + l) \JZ — 1 (n — a 2 ) p (2 p — 1) 
+ (?r — 2 p + a 2 ) n 


(32t t 2 Q 2 ) p T+^P-^b^ (Al n ~ 2 S 


^(cx*+2n-2p-l)/(n-l) y Z 

(a 2 + 1) y/Z=Tb a2 n f4l n ~ 2 S' ( "- q2)/( "- 1) 


(n-2p+a 2 )/(1—n)‘ 


1/(2p—1) 


^-2 v ^ (a2 - 2 " +1)/(n “ 1) 


16 7T / J 


= 0. 


(44) 


Considering S , Q and J as a complete set of extensive quantities for the mass M(S, Q, J), we should define conjugate 
intensive quantities to them. These quantities are temperature, angular velocities and electric potential 


T = 


/dM 

\dS 


, a = 


J ,Q 


/dM 

\dl1 


,U = 


S,Q 


fdM\ 

V J s, j 


(45) 


One can check numerically that quantities defined by Eq. (1451) coincide with Eqs. (1251) . (l26l) and (1421) provided C is 
chosen as C = (n — l)p 2 /n. Note that in the case of linear Maxwell field (p = 1), C reduces to 1 as we expect |27[ . 
Therefore one can conclude that thermodynamics first law 


k 

dM = TdS + Y^^idJi + UdQ, (46) 

2—1 

is satisfied. As one can see from U is proportional to value of A M at infinity and therefore U diverges if A^ 
diverges at infinity. On the other hand, it is obvious from thermodynamics first law that M is dependent on the value 
of U . Therefore, as we mentioned in section m, A M should be finite at infinity in order to have a finite mass. We 
took this fact into account for finding constraints on p and a (see Eq. ®). 


IV. STABILITY IN CANONICAL AND GRAND-CANONICAL ENSEMBLES 

In this section, we are going to study thermal stability of rotating black branes in the presence of a power-law 
Maxwell field in canonical and grand-canonical ensembles. It is necessary for a thermodynamic system to discuss 
its thermal stability. Thermal stability is investigated to ensure that the entropy of system is at local maximum or 
equivalently the internal energy of system is at local minimum. Therefore, the stability of a rotating black brane as 
a thermodynamic system can be studied in terms of entropy S(M, Q, J) or its Legendre transformation M(S,Q, J). 


















9 



FIG. 2: The behavior of (d 2 M/8S 2 )q, j versus a < 1 with l =si b = 1, q = 0.8, S = 1.25, n = 5 and p = 2. 



FIG. 3: The behavior of H^qj versus a < 1 with l = b = 1, q = 0.8, S = 1.25, n = 5 and p = 2. Note that curve corresponding 
to r+ = 1 has been devided by 10. 


In terms of mass M(S,Q, J), the local stability in any ensemble implies that M(S,Q,J) be a convex function of 
its extensive variables. Typically, this behavior is studied by calculating the determinant of the Hessian matrix of 
M(S,Q, J) with respect to its extensive variables X t , = [d 2 M/dXidXj] [|l|, (42j. H X(X > 0 guarantees 

that the system is thermally stable. The number of thermodynamic variables is ensemble-dependent. In canonical 
ensemble, the charge and angular momenta are fixed parameters and consequently the determinant of Hessian matrix 
H XiX . reduces to (d 2 M/8 S 2 )qj . Therefore, in order to find the ranges where the system is at thermal stability, 
it is sufficient to find the ranges of positivity of (d 2 M / 3S 2 )q^j where the temperature T is positive as well. In 
grand-canonical ensemble Q and J are not fixed parameters. 

We discuss thermal stability for uncharged and charged cases, separately. First we discuss the uncharged case. It is 
notable to mention that in the case of uncharged rotating black branes (d 2 M/dS 2 )j is exactly one which is obtained 
in H3: 


d 2 M \ = n ( a 2 + l) [(5 2 - l)(n + 1 - 2a 2 ) + 5 2 (l ^ a 2 )] (2 _ n _ Q 2 )/(a 2 +1) 
dS 2 J 3 7 rS 2 H-™&(™-3)'r [{a 2 + n - 2)S 2 + 1 - a 2 ] T+ 


(47) 


Therefore, in canonical ensemble we just briefly review the results of them. Since > 1, (1471) is positive provided 
a < 1. Therefore uncharged rotating black branes are stable in canonical ensemble for a < 1. Note that for uncharged 
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FIG. 4: The behavior of T versus a < 1 with l = b = 1, q = 0.8, S = 1.25, n = 5 and p = 2. 



FIG. 5: The behavior of {d 2 M / 8 S 2 )q ,3 versus a < 1 with l = b = 1, r+ = 1.5, E = 1.25, n = 5 and p = 2. 


case temperature is always positive (see Eq. (E5l) ~). In grand-canonical ensemble H|j 


can be calculated as: 


H 


M 

S3 ~ 


16 (1- + 
52(ra-l)7^4 ^ a 2 _|_ n _ 2) 2 2 _|_ 1 _ q.2J 


(48) 


From Eq. (l48l) it is obvious that > 0 provided a < 1, which is similar to the case of the canonical ensemble. 
From the above arguments we conclude that the uncharged rotating black branes are thermally stable provided a < 1 
in both canonical and grand-canonical ensembles. 

For charged case, we discuss the stability for a < 1 and a > 1, separately. Since charge does not change the stable 
solutions to unstable ones [43}, for a < 1 we have thermally stable charged rotating black branes. This fact is shown 
in Figs. O [3l [5j [6l El and HI Figs. [2] and [3] show that for a < 1, the obtained solutions are always stable in canonical 
and grand-canonical ensembles for any value of r+, as we expect. Since T > 0 guarantees that we have black branes 
for our choices, one should choose q < q ex t- The behavior of temperature is depicted in Fig. [I] with same constants 
as Figs. E] and El Behaviors of (d 2 M/dS 2 )Q .j and H^qj with respect to a < 1 for different choices of charge q are 
plotted in Figs. Eland El These figures once again show that charge cannot change stable solutions to unstable ones 
and therefore we have stable charged rotating solutions for a < 1 as we had in uncharged case. Fig. [T]shows positivity 
of temperature T for solutions that their thermal stabilities have been depicted in Figs. El and El Figs. El and E] depict 
stability of solutions for a < 1 for different values of p in canonical and grand-canonical ensembles respectively. The 
behavior of temperature for latter case is illustrated in Fig. [TOl For a > 1, one can understand from Fig. [TT] that 
event horizon radius of stable black branes encounter an upper limit r +max in both canonical and grand-canonical 
ensembles. The value of this upper limit is greater in canonical ensemble than grand-canonical ensemble. The effect 

























11 



FIG. 6: The behavior of Hggj versus a < 1 with l = b = 1, r+ = 1.5, H = 1.25, n = 5 and p = 2. Note that curve corresponding 
to q = 0.3 has been devided by 10. 



FIG. 7: The behavior of T versus a < 1 with l = b = 1, r+ = 1.5, S = 1.25, n = 5 and p = 2. 

of a on stability of solutions in both canonical and grand-canonical ensembles are depicted in Fig. [12] There is again 
an ensemble dependent upper limit this time on a that for values greater than it solutions are no longer stable. The 
value of a max is smaller in the grand-canonical ensemble. In terms of 1/2 < p < nj 2, stability is shown in Fig. 1131 In 
contrast with r+ and a, there is a lower limit for p i.e. for p > p m ; n , we have stable solutions. p m ; n is again ensemble 
dependent as well as r+ max and a max . For nj 2 < p < n— 1 where a has a p-dependent lower limit, numerical analysis 
confirm the result of investigation for 1/2 < p < n/2 , i.e. there is again a lower limit p m ; n that for values lower 
than it solutions are unstable. Stability in terms of charge is depicted in Fig. [TT] In this case there is an ensemble 
dependent q m in in each of ensembles that for q greater than it black branes are stable. The value of g m ; n is greater in 
grand-canonical ensemble. 


V. CLOSING REMARKS 

In this paper, we constructed a new class of higher dimensional nonlinear charged rotating black brane solutions 
in Einstein-dilaton gravity with complete set of rotation parameters. The nonlinear electromagnetic source was 
considered in the form of the power-law Maxwell held which guarantees conformal invariance of the electromagnetic 
Lagrangian in arbitrary dimensions for specific choices of power. Due to the presence of the dilaton Held, our solutions 
are neither asymptotically flat nor (A)dS. We showed that in case of power-law Maxwell Held, one needs two Liouville 
type potentials in order to have a rotating black brane solutions while in the case of linear Maxwell source just one 
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FIG. 8: The behavior of ( d 2 M/dS 2 )Q,j versus a < 1 with l = b = 1 , r+ = 1.5, E = 1.25, n = 5 and q = 0.8. Note that curve 
corresponding to p = 0.75 has been rescaled by the factor 1.5. 



FIG. 9: The behavior of Hggj versus a < 1 with l => b = 1 , r+ = 1.5, S = 1.25, n = 5 and q = 0.8. Note that curves 
corresponding to p = 0.75 and p = 2 have been rescaled by the factors 10 and 10 _1 respectively. 


term is needed [ 27 |. The extra dilaton potential term disappears for p = 1 . All our results reproduce the results of 
| 27 ( in the case of linear charged rotating solutions where p = 1 . 

Demanding from one side that the value of total mass should be finite and from the other side that the effect of 
mass term in metric function should disappear at infinity, we found some restrictions on p and a. The allowed ranges 
of these two parameters is as follows. For 1/2 < p < n/ 2 , we have 0 < a 2 < n — 2 while for n/2 < p < n — 1 , we have 
2p — n < or < n — 2. For these permitted ranges, our solutions are always well-defined. Also, in these ranges of p 
and a, the charge term in metric function /(r) is always positive and dominant in the vicinity of r = 0. Therefore, 
Schwarzschild-like solutions are ruled out. However, solutions with two inner and outer horizons, extreme solutions 
and naked singularities are allowed. 

In order to study thermodynamics of charged rotating black branes, we calculated mass, charge, temperature, 
entropy, electric potential energy and angular momentum. Using these quantities we obtained Smarr-type formula 
for the mass M(S,Q, J) and showed that the first law of thermodynamics is satisfied. Next, we analysed thermal 
stability of the solutions in both canonical and grand-canonical ensembles. These investigations showed that for a < 1 
charged rotating black brane solutions are always stable with any value of the charge parameter. For a > 1, there is 
a r+max in each of ensembles that we have stable solutions provided their radius are smaller than r+ max . In terms 
of a(> 1), the solutions is changed from stable ones to unstable ones when they meet an a max . Versus p and q 1 we 
showed that the solutions encounter p m i n and q m i n respectively so that solutions with q and p parameters lower than 
them are unstable for a > 1. All r +max , a max , q m i n and p m i n values depend on the ensemble. 
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FIG. 10: The behavior of T versus a < 1 with l = b = 1, r+ = 1.5, S = 1.25, n = 5 and q = 0.8. Note that curve corresponding 
to p = 0.75 has been rescaled by the factor 1.5. 



FIG. 11: The behavior of T (solid curve), 10(<9 1 2 AF/<9S' 2 )q,j (dashed curve) and 10 3 4 Hm, (dashdot curve) versus r+ with 
l = b=l, , = 0.5, a = 1.35, S = 1.25, n = 5 and p = 2. 


It is worth noting that the higher dimensional charged rotating solutions obtained here have flat horizon. One may 
interested in studying the rotating solutions with curved horizon. Specially the case of spherical horizon will be a 
good extension of Kerr-Newmann solution. It seems that the study of the general case is a difficult problem. However, 
it is possible to seek for slowly rotating nonlinear charged solutions with curved horizon. The latter is in progress. 
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